DEFORMATIONS OF THE MONGE /RIEM ANN HIERARCHY AND 
APPROXIMATELY INTEGRABLE SYSTEMS 
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Abstract. Dispersive deformations of the Monge equation u u = uu x are studied using ideas 
originating from topological quantum field theory and the deformation quantization programme. 
It is shown that, to a high-order, the symmetries of the Monge equation may also be appropri- 
ately deformed, and that, if they exist at all orders, they are uniquely determined by the original 
deformation. This leads to either a new class of integrable systems or to a rigorous notion of 
an approximate integrable system. Quasi-Miura transformations are also constructed for such 
deformed equations. 



1. Introduction 

Consider a general scalar evolution equation 

u t = K(u) (1.1) 

where K (u) is a smooth function depending on u and the z-derivatives of u, but not explicitly on 
the independent variable x . A characteristic feature of an integrable system is the existence of an 
infinite hierarchy of compatible, commuting flows, these being the generalized symmetries of the 
original equation ( [2.3| ). Such symmetries have the form 

u T = Q{u) 

(where again Q is a smooth function of u and the x-derivatives of it, but not explicitly on the 
independent variable x) and must satisfy the equation 

K'Q - Q'K = 

where the prime denotes the Frechet derivative 

K'{u)Q=^-K{u + eQ) 

The paradigm for such constructions is the Korteweg-deVries equation 

u t = uu x + u xxx . (1.2) 

For this equations, the hierarchy of symmetries may be constructed by exploiting the bi-Hamiltonian 
structure of the KdV equation. Given the Hamiltonians 

Hi = d , H 2 = d 3 + ^ud + ^u x 

one may form the recursion operator 1Z = 7i 2 'H^ 1 with which the symmetries take the form 

u Tn = K n u x (1.3) 

(for more details see, for example, H). Such flows constitute the KdV hierarchy. 

There has been much interest recently in the properties of the dispcrsionless limits of integrable 
systems, these being examples of equations of hydrodynamic type. As well as being interesting 
integrable systems in their own right, they play an important role in topological quantum field 
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theories and the construction of higher-genus Gromov-Witten invariants [DZ1, DZ2|. Under the 
scalings 

t i ^ fi _1 f 

x i— > fi a: 
U I— > u 

the KdV equation scales to 

ttf = uu>x + fc 2 u xxx 

and under the classical limit ft — > one obtains the Monge equation 

u t = uu x . (1.4) 

This singular limit has a drastic effect on the solution space of the KdV equation, for example 
soliton solutions to not survive, but many properties do survive, notably the form of conservation 
laws and the existence of a bi-Hamiltonian hierarchy. Thus a hierarchy of symmetries of the Monge 
equation may be obtained by scaling the KdV hierarchy and taking the h — > limit. This results 



in the Monge hierarchy (also called in [DZ2| the Riemann hierarchy) 



However, the Monge equation has many more symmetries; given any suitably differentible function 
f(u) then 

u r = f{u)u x (1.5) 

is a symmetry, labelled not by a discrete integer, but by a C 2 (R)-function. 

The purpose of this paper is to study not the dispersionless limit but the reverse, regarding 
the KdV equation as a deformation of the Monge e qua tion and constructing, term-by-term, the 
corresponding deformation of the general symmetry ( |l.5|) . In fact a slightly more general equation 
than the KdV equation will be studied, namely the equation 

u t = uu x + (a(u)u xxx + b(u)u xx u x + c(u)w|) 

where a , b , c are arbitrary functions. The motivation for this comes from the work of Eguchi et 
al. [EYY], Dubrovin and Zhang |DZ1, DZ2] and, in particular, Lorenzoni Q. Their approaches 



utilize the bi-Hamiltonian structure of their systems. No such assumption will be made here and, 
as a result, certain of Lorenzoni's results will appear as a special case of the constructions that 
will appear in the subsequent sections. 

2. Notation and the form of the general problem 

Let T> denote the linear vector space/ring of differential polynomials in \ui,U2, •••} where 
ui = u x ,U2 = u xx etc., with coefficients in the space C°°(K) of infinitely differentible functions, 

2J = C , ~(R)[tti J « a ,...], 
so an arbitrary element F ET> takes the form 

F = aj(u)uj 

\I\«x> 

where ai(u) £ C°°(M.) and / is a multi- index. To each monomial one may assign a scaling degree 
by counting the number of x-derivatives it contains, so 

deg(uk) = k , k > 1 , 
deg(AB) = deg(A) + deg(B) , A , B G V 

and no degree is assigned to the coefficient functions (this is very similar to the ideas used in the 



theory of normal forms [KM], but a different grading is used there). With this, V decomposes 



DEFORMATIONS AND APPROXIMATELY INTEGRABLE SYSTEMS 



3 



into a sum of sub-spaces/rings consisting of terms of the same degree 

oo 

2> = 0£> fc) . 

Thus, for example, 



fc=0 



T> 3 = span{u xxx , u xx u Xl u 3 x } . 

Clearly dirnVk = number of partitions of k , denoted by #(&)■ With this notation it is easy to 
describe the problem to be studied. 
Consider the evolution equation 

u t = Ki (2.1) 
where K\ — uu x S T>\ and corresponding symmetry 

u T = Qi (2.2) 



where Q\ — f(u)u x G T>\ . Suppose the equation (2.1) is deformed by an arbitrary element 
K 3 6 V 3 : 

u t = K[h], 



Can one construct a corresponding deformation of the symmetry (2.2) 

u T = Q[h] , 



— /J h 2n Q2n+l 



n=0 

with Qi 6 T>i and Q\ as above? The deformation parameter h should be regarded as a formal 
parameter, labelling the scaling dimensions of the terms it precedes, rather than a small parameter. 
One may, at the end of all calculations, set h = 1 . 

By definition, the symmetry Q[h] must satisfy the governing equation 

K[h)'Q[h] - Q[h]'K[h] = o 

and equating powers of h gives the basic equation 



K[Q 2n +l ~ Q'2n+1 K 1 = Q2n-1^3 ~ K' 3 Q 2n -l 



neN. 



(2.3) 



By construction the equation for the lowest power h° is automatically satisfied. This equation (2.3) 
looks complicated but, since all the elements lie in various spaces with fixed scaling dimension, it 
reduces to an over-determined linear system. 

3. INTEGRABLE DEFORMATIONS OF THE MONGE EQUATION 



In this section equation (2.3) will be studied in more detail, with 
Ki = uu x , 

K 3 = a(u)u xxx + b(u)u xx u x + c(u)u x , a,b,c£ C c 
Qi = f{u)u x , 



so 



u t = uu x + h 2 (a(u)u xxx + b(u)u xx u x + c(u)u^) . 
With these (2.3) becomes 

(a'u3 + b'u 2 ui + c'u\)+ 



'2n+l 



d 



^ [(uux) r - uu r+ i]-g— - ui 



Q271+1 — 



(ad 3 + b Ul d 2 + bu 2 d + Zu\d)- 



12n-l 



(3.1) 



(3.2) 



4 



I.A.B. STRACHAN 



or schematically as 

n(Q 2 „+i) = S(Q 2n -i) (3.3) 

where 

S : V n ► £>n+3 ■ 

One may regard this equation in a number of ways: 

• as a first order quasi-linear equation in the independent variables u\,u 2 etc., and solve it 
using the method of characteristics. However, one requires, for the purposes of this paper, 
the solution to lie in T>2 n +i which is hard to guarantee using this method; 

• as a problem in differential Galois theory - one requires the solution to this linear equation 
with coefficients in a particular ring to lie in that ring and not in some extension; 

• as an overdetermined linear system. 

This third approach is the one that will be using here. This interpretation holds because d u does 
not appear in the operator 7Z, so no derivatives of the unknown coefficients appear. By introducing 
a basis for T>2 n +3 one obtains #(2n + 3) linear equations for #(2n + 1) unknowns. 

In the next section it will be shown that the rank of 1Z (viewed as a matrix) is maximal, and 
hence the solution to this linear problem, if it exits, is unique. Existence is more problematic. 
While the form of the matrix TZ is easy to understand, the 'source' vector S(Q2n-i) is complicated, 
depending on all lower order solutions, and so it would be extremely difficult to compute, in general, 
the rank of the extended matrix. By direct calculation (using Mathematica) explicit solutions up 
to 0(H 6 ) have been calculated. This suggest the conjecture that a solution exists for all n, so a 



formal symmetry of equation (3.1) of the form fi 2n Q2n+i & 2? exists. 

4. Uniqueness 



One may to prove that the equation (3.3) has a unique solution is to introduce an ordered basis 



for each subspace T>k of monomials uj - the reverse lexiographic ordering - so, for example 

V 5 = span{u 5 , UiUi, u 3 u 2 , u 3 ul, u 2 U\,U2U 1 , u{\ 
and this basis may be order 

M5 >- u^ui >~ U3U2 >- u 3 u\ >- u\u\ >- U2u\ >- u\ . 

The symbol max{v} will denote the highest basis vector in the expansion of v in this ordered 
basis. 

Lemma 1. In the reverse lexiographic ordered basis the operator 7Z : T> n ► Z?„+i is lower- 
triangular. 

It then follows that the equation IZv = has the unique solution v — and hence that the solution 



to (3.3), if it exists, is unique. 



Proof Let ui be a monomial in T> n , / being a multi-index. From the explicit form of 1Z it follows 
that 

max{7Zui} = uiU\ . 

Suppose that uj >- uj . Then 

max{TZui} — uiu\ >~ uju\ = max{TZuj} . 
Hence, in this basis, the matrix 1Z is lower-triangular. 

For example, 

TZ[giu xxx + g 2 u xx u x + g 3 ul] = 3giu xxx u x + 3giu 2 xx + 3g 2 u xx u 2 x + 2g 3 u x 
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so in the reverse lexiographic basis, 1Z may be represented as the 3x5 matrix 



which clearly is of maximal rank. 



K 



/0 0\ 

3 

3 

3 

\0 2 J 



5. Existence 

Recall that the given data is K\ , K% and Q\ . Writing 

Qz = 9\u xxx + g2u xx u x 

and solving the over-determined linear system yields 

.9i = af , 

g 2 = bf' + 2af", 

93 



93U X 



cf + \bf" + \af" 



Thus a solution exists. This gives an approximate symmetry, valid to 0(h 2 ) . Recycling this 
solution gives 



Qs = hi ub + h 2 U4U1 + h 3 U3U2 



U3U1 + h$ u\u\ + ha U2u\ + 



h 7 u\ 



(5.1) 



where the hi are explicit functions of a, 6, c and / and their derivatives. These are given in the 
appendix. The 0(h 6 ) have also been calculated, again using Mathematica, but little purpose is 
served by presenting them here - it suffices to say that they exist. Thus for arbitrary functions 
a,b,c a unique symmetry, labelled by a suitably differcntablc function / exists up to 0(h 6 ) . The 
existence of solutions up to this order of the over-determined linear system suggests the following: 



Conjecture 1. A unique formal symmetry to the generalized KdV equation (3.1) exists, labelled 
by a C°°(R)- function. 

If the conjecture is false, then this would raise the questions: 

Conjecture 1'. At what order does the above conjecture fail? Is there a way to determine this 



order 



a priori 



i? 



Note that if the conjecture is true then one obtains an entirely new integrable hierarchy, depending 
on three arbitrary functions and with flows labelled by another arbitrary function. If it is false, 
then one may obtains the notion of an approximately integrable system, which has an infinite 
number of conservation laws up to some fixed order. Either outcome is of interest. In a particular 
case, Lorenzoni jO has numerically observed elastic soliton scattering indicative of integrability. It 
would be interesting to see if basic integrability results such as inverse scattering could be modified 
(by, say, only including terms up to a given order) to include such approximate integrable systems. 

5.1. Formal symmetries of the KdV equation. In this approach, the KdV equation itself is 
recovered as the special case a = 1, b = c = , and the symmetry takes the form 



f(u)u x 



f'u xx + \ulf" 



lf {5) < + l ^uiu xx + A/(^L + If 



10 



10" 



(3) u u + _f(2) u 

^a; Ul XXX 1 ~ .1 Ul xxx:i 
5 



+ 0(k 6 ). 



(again the 0(h 6 ) terms have been calculated explicitly, though not displayed here). One may 
prove from the now simplifier version of (pj) that for the series to terminate the function / must 
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be a polynomial, say / 



r, a r u r , and one recovers a linear combination of the flows obtained 



via the recursion operator fll.g), 



Vr=0 



It would be of interest to see if the general flow generated from an arbitrary / could be written 
as u T — F(lZ)u x for some suitable F . 

Much has been written about the symmetries of the KdV equation. However in such approaches 
the order of the symmetry is fixed, and the lower-order terms are determined from the higher ones. 
Here is approach is opposite - determine the higher-order terms from lower ones. The symmetries 
obtained in this was will turn out to be formal - infinite series. This approach is motivated, as 
was mentioned in the introduction, by ideas originating in topological quantum field theories, 
where the deformation is known as a genus expansion, and from the deformation quantization 
programme, where in both areas one constructs higher-order terms from lower-order ones. 

5.2. Local Hamiltonian systems. An important subclass of systems in this class are those that 
are Hamiltonian with respect to the local Hamiltonian operator Hi = d , so the generalizes KdV 
equation (fO) takes the from 



Ut 



ou 



for some Hamiltonian H . This places the following restriction on the arbitrary functions a, 6, c : 

a (it) = s(it) , 
= 2s'(u), 



b(u) 
c(u) 

where s is an arbitrary function. With these 
Ut = uu x + h 2 



\s"{u) 

ll) may be written in Hamiltonian form 



f 



s{u)u xxx + 2s'(u)u xx u x + ^s"(u)ul 



d_6_ 

dx Su 



u 

~6~ 



The formal symmetries arc also Hamiltonian: 



d_S_ 

dx Su 



^ { ±s 2 fu 2 xx i 



0(h e ) 



10" •' itx 4! 

(where 9 2 /'~ 2 ) = /). Again, the 0(ft 6 )-terms have been calculated, and are also Hamiltonian, 
which suggests the following: 

Conjecture 2. The formal system in Conjecture 1 is Hamiltonian to all orders. 

This system has been extensively studied, for /(it) = it", by Lorenzoni [0. He showed that 
the systems is, up to 0(fi 4 ), bi-Hamiltonian, with the second Hamiltonian structure being a 
deformation of the second Hamiltonian structure of the Monge/Riemann hierarchy. The existence 
of such terms in this deformation is controlled by a certain cohomology group. 



6. Trivial and quasi-trivial Miura transformations 



In | DZ2 1 transformation of the form, 



J2h k F k ( 



u;u x , 



,u 



k=l 



where are quasi-homogeneous rational functions, were considered together with the correspond- 
ing action on bi-Hamiltonian pencils. Such transformations were called quasi-Miura transforma- 
tions. In particular the notion of trivial and quasi-trivial transformations were given. In the 
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context of this paper, where local bi-Hamiltonian structures are not consider, an evolution equa- 
tion u t = uu x + h r K r [u] will be said to be quasi-trivial if it transforms under a quasi-Miura 
transformation to the Monge equation vt = vv x , and trivial if the further condition that the 
functions Fk are polynomial is satisfied. 
With the ansatz 



4 r -2 F r [u,Ui, . . . ,U 6r - 2 ] , F r €l>6r-2 

one may easily obtain recursion relations for the F r which take the form 

lZ(F n ) — function of F\ , . . . , 77 n _i 

where 

7c : T> n — ► 

is given by 

K = uiK+ (2 - An)u\ 

with 7c being the previously introduced recursion operator. It follows immediately from a repeti- 
tion of the argument in section 4 this that the quasi-Miura transformation, if it exists, is unique, 
since the operator 7c is lower triangular in the reverse lexiographic ordered basis. 



h 2r 



Lemma 2. The generalised KdV equation (3.1) is, up to 0(h ), quasi-trivial. 



This may be proved by direct calculation. The quasi-Miura transformation being 

h 2 (\ 1 \ h 4 

v = u + — I -a(u)(u xxx u x - u 2 xx ) + -b(u)u xx u 2 + c(u)u x ) + — f(u, u x , . . . , u Wx ) 



where / G 2?io has been explicitly calculated, though not displayed here. In [DZ2] it was shown 
that the KdV equation is quasi-trivial to all orders, this proof using the existing machinery for 
the KdV equation. Conjecturally the generalized KdV equation studied here will be quasi-trivial 
to all orders, though to prove it would involve the development of a lot of associated machinery. 



One immediate question is whether the quasi-trivial equation (|3.1| ) can, for suitable choice of 
functions a , b , c be trivial. From the above form of the quasi-Miura transformation it follows that 
at 0(h 2 ) one requires a(u) = and at 0(h 4 ) one requires, in addition, b(u) — 0. After this no 
new constraints appear, at least to 0(7i 10 ) , the explicit quasi-Miura map being 

v — u — h 2 cu 2 + h 4 ^—cc' uf + 2 c 2 ufu2 

^ ( -28cc' 2 ui 6 14c 2 c"wi 6 2Qc 2 c'u 1 4 u 2 , q , 2 4c 3 ui 3 w 3 ' 

+h b — - - Ac 6 ui u 2 — 1 

^15 15 3 3 

tR /836cc' 3 wi 8 1268 c 2 c'c" ui 8 16c 3 c ( 3 )ui 8 1138 c 2 c' 2 m 6 u 2 

+h s — H — H — H — - 

\ 315 315 35 45 

128 c 3 c" Ul e U2 112 c 3 c'u 1 4 u 2 2 



15 

o 4 4 
2 C U\ Ui 



8 c 4 ux 2 u 2 3 + 8 c 3 d iti 5 !i 3 + 8c 4 Mi 3 u 2 u 3 + 



+ o{h w ) 



3 

Thus the first order, third degree evolution equation 

u t = uu x + h 2 c(u)u x 

is trivial to 0(h s ) and conjecturally to all orders. Further properties of this system are outlined 
in the next section. 
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7. First order deformation of arbitrary degree 
Exact results to all orders may be obtained in the above special case a = b = 0, c(u) = a(u) . 



The basic equation (3.1) simplifies to 

u t = uu x + h 2 a{u)ul , aeC°°(R) (7.1) 
and the governing equation (|2.3| ) simplifies drastically. With the ansatz 

one obtains the recursion relation 

n/3 n = (1 - n)c//3 n -i + afi' n _ x , n>l, 
and hence, since (3 = (3o is arbitrary, the solution 



Thus the general symmetry of the first order, third degree equation (7.1) is 

n! 



Thus the space of commuting flows for equation (7.1) is labelled by an arbitrary function (3 6 
C°°(M). This series truncates if and only if j3 is a polynomial. This symmetry is, within the 
constraints of this paper, unique. It can easily be checked that all the symmetries commute 
amongst themselves, i.e. 

for arbitrary C°°-functions a, /3, 7 . 



The equation (7.1) is Hamiltonian with respect to the non-local Hamiltonian operator 

H = u x d~ x u x , 

so 

u t = (li.— ) f ( u ~ h 2a ( u ) u x ) dx . 



' Su 

The symmetry itself is also Hamiltonian 



^ = H ~Su~ 



where the Hamiltonian is 

H„ = 



r=0 ^ ' 



and these are all conserved with respect to all the symmetries, as one would expect for an integrable 
system. Owing to the arbitrariness in the functions a and /3 is unlikely that ([O]) will possess the 
Painleve property. 

8. Conclusions 

The central philosophy behind this paper is that the dispersionless limits of integrable systems 
should not been seen as some limit, but as the kernel of the dispersive hierarchy itself, encapsu- 
lating its integrability. To describe the full dispersive hierarchy one just has to specify a suitable 
'deformation vector' in some space of possible first-order deformations, from which the dispersive 
hierarchy, if it exists, may be uniquely reconstructed. In this paper this vector is K3, living in the 
space of possible first-order deformations P3 . If the above conjecture is true, then this specifies 
the hierarchy at all orders. For equations of hydrodynamic type constructed from semi-simple 



Frobenius manifolds this choice is automatic [DZ1, DZ2| - elliptic Gromov-Witten invariants may 



be constructed from the genus-zero data. See also [ EYY] , where examples are given where higher- 
order deformations do not exist. 
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As an multi-dimensional example consider the system 

u t = U o u x , 

where o is the product is some A-dimensional, commutative A-dimensional algebra. With ceratin 
simple additional conditions on this algebra one may show that a hierarchy exists if and only if 
the algebra is a Jordan algebra Q . It then may be shown that the deformed system 

u ( = Uoua h 2 u xxx 

(corresponding to the deformation vector 

Uxxx e T>£ = span{u xxx , u xx u J x , u l x v? x u k x , k = 1, . . . , N} 

where U = u l e,i and the form a basis for the algebra) defines a full dispersive hierarchy with no 
extra conditions required. It would be interesting to repeat these calculations for a more general 
vector in - the results in this paper being for N = 1 only. 

It may also be of interest t o ge neralize the construction, in particular the quasi-triviality prop- 
erty, to include ifs-terms in ( |2.3| ), 

u t = Ki + h 2 K 3 + h 4 K 5 . 

A number of known systems fall into this class (with K\ — u 2 u x ), including the higher-flow of the 
KdV equation, KdV$ , the Sawada/Kotera equation, the Caudrey/Dodd/Gibbons/Kaup equation 
and the Kaup/Kupershmidt equation. This paper suggests a possible classification of such fifth- 
order systems. Also, little mention has been made of bi-Hamiltonian structures in this paper. 
Such structures may be generated from the known bi-Hamiltonian pencil of the Monge equation 

{u(x),u(y)}\ = (u- X)S'(x -y) + -u x S(x - y) 
via the appropriate quasi-Miura transformation. These, however, are unlikely to truncate. 

The results outlined here have been computational; to obtain rigorous results at all orders in 
the expansion new techniques will have to be developed. In particular, a better understanding of 
the over-determined linear systems from which results of this paper were obtained will be required. 

Appendix . The 0(h 4 ) order terms 
The explicit form of the hi in equation ( |5.l|) are: 

h = (3a 2 /")/5, 

h 2 = a6/" + (8aa7")/5 + (9a 2 / (3) )/5, 
h A = 2abf" + 2aa'f" + 3a 2 / (3) , 

h 4 = (2b 2 f")/5 + {llacf")/5 + ba'f" + {13ab'f")/W + (7af"a")/5 + 

(5a6/ (3) )/2 + (37aa'/ (3) )/10 + (23a 2 / (4) )/10 , 
h 5 = (8b 2 f" + 5b(a'f" + 7a/ (3) ) + a(Ucf" + 26b' f" + Ufa" + 44a' '/^ + 31a/ (4) ))/10 , 
h 6 = (% 2 /^+&(18c/" + 967" + 4/"a" + 15a7 (3) +25a/( 4 )) + 

2{Aca'f" + Wac'f" + 7af"b" + I6acf^ + Uab'f^ + 12aa"f {3) + 

2af"a^ + 18aa7 (4) + 8a 2 /< 5 )))/10 , 
h 7 = (4c 2 f" + b 2 f^ +2c(b'f" + 2bf i3) +a'/ (3) +2a/ (4) ) + 

b(2c'f" + f"b" + 2b' + a"/ (3) + 2a' + 2a f&) + 

a(2f"c" + 4c' + 3b" + f^aW + f'bW + 3&7«> + 3a" + 3a'/ (5) + af^))/8 
The fifteen coefficients in the next-order term have been calculated, but are not presented here. 
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